Anomalous features of the proximity effect in triplet superconductors 
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Anomalous features of the proximity effect specific to diffusive normal metal / triplet superconduc- 
tor (DN/TS) junctions are studied within the quasiclassical Green's function formalism. The pair 
amplitude /jv as a function of energy e in DN, satisfies an anomalous relation, /jv(e) = — /jv( — e), 
contrary to that in singlet superconductor junctions case where /jv(e) = fpf{~ e )- Such an unusual 
e dependence is responsible for a zero energy peak in the local density of states and is a source of an 
anomalous penetration of the applied magnetic field into DN. Our results are relevant to Sr2Ru04, 
which is considered to have sl p x + ip^-wave symmetry. 



Nowadays spin-triplet superconductivity is a fascinat- 
ing and an important issue in condensed matter physics. 
The existence of triplet superconductivity has become 
promising after a series of experiments in S^RuO^^. 
Other possible triplet superconductors are UPts 3 and 
(TMTSF^X 4 . Phase sensitive phenomena— in junctions 
involving triplet superconductor is an important and rel- 
atively unexplored subject. One open question is how 
the proximity effect manifests itself in junctions between 
triplet superconductors and normal metals. 

The proximity effect can be described in terms of the 
coherence between electrons and Andreev reflected holes 
in a normal metal. In unconventional superconductors, 
the internal phase reflecting sign change of a pair poten- 
tial significantly affects charge transport since this sign 
change is a source of a midgap Andreev resonant state 
(MARS) at the surface of a superconductor. The origin 
of the MARS& is the interference effect of a quasiparti- 
cle at the interface^. In ballistic case, MARS leads to 
the formation of zero energy peak (ZEP) in local density 
of state (LDOS). However, the interference is destructive 
for (i-wave superconductor /diffusive normal metal (DN) 
junctions^ and therefore no ZEP is expected in LDOS 
in DN. On the other hand, in diffusive normal metal 
/ triplet superconductor (DN/TS) junctions the MARS 
penetrates into DNltt and ZEP occurs in the LDOS in 
DN. 

Recently, a search of new phenomena in triplet super- 
conductors has attracted a lot of attention. For example, 
recent experiments indicated possible existence of MARS 
in the tunneling spectroscopy of Sr 2 Ru04 and Ru micro 
inclusion system 1 ^. Therefore, it is important to explore 
the physics of proximity effect in DN/TS junctions and 
to address various situations relevant for the actual junc- 
tions. 

In the present paper, we study the anomalous proxim- 
ity effect in N/DN/TS and TS/DN/TS junctions (N de- 
notes normal electrode) in the framework of the Keldysh- 
Nambu quasiclassical Green's function formalism estab- 
lished in previous work«2i±fliii. We find a crucial rela- 
tion in pair amplitudes in DN, /jv(s') = — /jv( — e ) m 



contrast to a relation /jv(e) = /at( — e ); which holds in 
DN attached to a singlet superconductor. Here e is the 
quasiparticle energy measured from the Fermi level. The 
function /jv(^) is commonly used in the quasiclassical 
Green function theory and characterizes the pair ampli- 
tudes in DN due to the proximity effects 1 ^. The property 
/jv(e) = —f^(—e) is responsible for the penetration of 
MARS into N/DN/TS junctions which provides the ZEP 
in LDOS. Moreover, we predict an unusual response of 
DN to an external magnetic field. We show that the 
magnetic field spatially oscillates in DN instead of usual 
Meissner screening. These unusual effects do not exist 
in the conventional proximity effect in singlet supercon- 
ductors. We also show that in TS/DN/TS junctions, the 
proximity effect is suppressed when the external phase 
difference across the junction tp is due to the destruc- 
tive interference of MARS's from the two superconduc- 
tors. On the other hand, at ip — n the constructive in- 
terference bridges the two MARS's in DN. 



The model of the system is shown in Fig. 1. We 
consider a junction consisting of a normal metal and 
a superconducting reservoir connected by a quasi-one- 
dimensional diffusive conductor (DN) with a length L 
much larger than the mean free path. The interface be- 
tween DN and TS has a resistance Rb while the DN/N 
interface has a resistance Rb> . We assume flat inter- 
faces and the insulating barriers at the two interfaces 
are modeled by the delta-functionii. We start from the 
case where the pair potential is chosen to be p x + ip y - 
wave pairing symmetrjii^ which is the promising candi- 
date of superconductivity in S^RuO^. As regards the 
spin structure of the Cooper pair, we choose S z = . 
The positions of the DN/N interface and the DN/TS in- 
terface are denoted as x = and x — L, respectively. 
We calculate the Keldysh-Nambu quasiclassical Green's 
functions by using the general boundary conditions given 
in Eqs. (4) and (5) in Ref.[lll To calculate the LDOS in 
DN, we focus on the retarded part of the Nambu-Keldysh 
Green's functions. For p x + ip y -w&ve pairing symmetry, 



2 



the quasiclassical Green's function in TS is given by 



9T~3 



(1) 



Rl± = fl±T\ + $2±T2 

with f 1± = Ai,±/ a/Aq - e 2 , / a± = -A 2 ,±/VA 2 -£ 2 , 
<7 = e/ \/e 2 — A 2 , where Ai± and A 2 ± are given by 
Ai± = ± A cos <f> and A 2 ± = Aosin^, respectively. Here 
Ai,+ + iAa,+ and Ai_ 4- ?A2 ; _ are the effective pair po- 
tentials for a quasiparticle with an injection angle (j) and 
(tt — 4>)i respectively. The Green's function in DN is ex- 
pressed by 



Rn{x) — sin0f 2 + cos0T3 



(2) 



with the proximity parameter 9. Since the T3 component 
of the matrix current in Eq. (4) of Ref. ifllh vanishes due 
to the absence of the Josephson current after the angular 
average of the matrix current as a function of ip, the fx 
component of the Green's function in DN is vanishing. 
As discussed in Ref. iflll) . function 9 in DN is determined 
by the Usadel equation 
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D——9 + 2iesin9 = 
ox 1 

supplemented by the boundary conditions 
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where 
Fx = 
F 2 = 
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dx 



(Fx), 
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J-2^2 f2 
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2 - T' + V cos 9 







(3) 

(F2), 
(4) 

■, (5) 
(6) 



with A = (cos^l/ 2 — sin^Lff), B = + cos 9^g), 

T x = T /(2 - + 2^l~%) and T' = T£. 

Here D is the diffusion constant i DN and (...) denotes 
the average over angle <f> defined as follows 

/7T/2 /•7T/2 
d<j> cos 4>Fx(2){(j>) I \ dcjyT^ cos 0, 

(7) 

is the resistance of DN, 0£ and #0 are the values of 
at x = L and x = 0, respectively. In the above, /1 = 
/i+ = —fx- and /2 = fi+ = /2- are satisfied. The 
transparencies at the DN/TS and N/DN interface are 
given by 



T A = 



4 cos 2 



Z 2 + 4( 



T' = 



4 cos 2 



Z' 2 + 4cos 2 </>' 



(8) 



where Z and Z' denote the barrier parameter at the in- 
terfaces. 

At first we discuss the unusual e dependence of 9. It 
is easy to confirm relations A( 2 )(— £,</>) = fx(2)( E i 
<?(-£, 0) = flr*(e,0), fi{e,-(f>) = fx(s,-(j>), h(e,~4>) = 



— /2(e, — 0), and g(e,—(f>) = g(e,<j>) by the definition of 
/1.2 and g. By applying these relations to Eqs. Q-lJSJ, 
we immediately find 



?*(£). 



(9) 



It can be shown from Eq. (22) in Ref. Q that this equa- 
tion holds for triplet junctions when time reversal sym- 
metry is preserved—. We note that 9(~e) — 9*(e)^ is 
satisfied in any singlet junction which preserves the time 
reversal symmetry. Thus, Eq. is an important and a 
unique relation in triplet junctions. This relation leads 
to unusual properties of the proximity effect in triplet 
junctions. In what follows, we explicitly describe x in 
p(e,x) and f^(e,x) — sm9(e,x). In Fig. 2, the local 
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FIG. 1: Schematic illustration of the model of the N/DN/TS 
junction with p x + ip y -wa,ve pairing 

density of states p(e,x) = Real[cos#(£, x)} and fw(e,x) 
with 9(e, x) = 9(e) = 9 at the center of DN (i.e., x = L/2 
) in N/DN/TS junctions arc plotted for p x and p x + ip y , 
where Z = 1, Z' = 1, and we fix a ratio Rd/ Rb at 1. The 
interface resistance at N/DN (Rb') is a parameters of the 
calculation. In order to achieve the convergence of LDOS 
in numerical calculation, we replace £ with e + ij where 
7 is small parameter, which we choose to be 0.005Ao. 
As shown in Fig. 2(a), the LDOS has a ZEP, which indi- 
cates that MARS due to the proximity effect exists deep 
in the DN. The height of the zero energy peak (ZEP) 
of LDOS increases with the decrease of Rd/Rb'- This 
is because the pair amplitude fj^(s,x) is more strongly 
confined within DN for larger Ry. In p^-wave pairing 
symmetry, the same tendency is found in p(e,L/2) as 
shown in Fig. 2(b). In s-wave junctions, in contrast to 
the triplet case, the amplitude of LDOS around e = 
decreases with the decrease of Rd/Rb' as shown in (c). 
For small Rd/Rb' , quasiparticles in N are separated from 
those in DN and the minigap exists in DN. With the in- 
crease of Rd/Rb' the states within the minigap arc filled 
in and the LDOS around e — increases. 

It is possible to understand the unusual e dependence 
of LDOS for N/DN/TS junctions with p^-wave symmetry 
in the limit Rd/Rb' —* 0. We concentrate on the case with 
< £ << A and << I. The spatial dependence of 
9 has a form 9 = 9 L0 + (x - L)6x/L + (x - L) 2 9 2 /(2L 2 ). 
Substituting fx = A cos 4>/ \/A 2 cos 2 (p — e 2 , fi = 0, and 
g — e/ ^Je 2 — A 2 , cos 2 4> with e — e + ij into Eq. JSJ, we 
find < Fx >^ Ao(sin#L + «cos#lo)/7- Further, using 
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Eqs. iJSJ and @, we find cos#lo and p(s,x 
cos 8 L0 



A R d D 



p{e,x) 



A R d D 



(10) 
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It is remarkable that the a;) is proportional to 
-iJl/Jsj). This unusual e dependence is specific to 
DN/TS junctions. 



- 1 (a)/?.v+/pv ' 


: 


(b)p-. 1 


a 


3 








a/ \ 


6 
















w 














Q. 

2 














i h- 


• 


£.-sA' 









\"B 7 
V 7 a 








1 

E/£ra 




1 


1 




s 


^ 




rvT i 
















V b _ 










e/En 



1 a-i 

5 



o 



FIG. 2: In the upper panel, local density of states p(e,L/2) 
is pfotted for N/DN/S junctions for various superconduct- 
ing state with Z = 1, Rd/Rb = 1, Z' = 1 and Z?Th = 
0.25Aq, respectively. (a)pa; + ip y -wa,ve, (b)p x -wa.ve and (c)s- 
wave pairing, a, Rd/Ry = 0.01; b, Rd/Ry = 1; and c, 
Rd/Rb' — 100. The corresponding pair amplitude /jv(e,L/2) 
with Rd/Ry = 1 are plotted for (a)^ + ip y -wave, (b)p x -wa,ve 
and (c)s-wave pairing, respectively, a: R,eal[/jv(e, £/2)], b: 
Imag[/ JV (£,L/2)]. 



As a result, the magnetic field in the DN behaves as 
H(x) ~ exp(— x/X(x)) with the local penetration depth 
A(x), which is given by 



1 TV sin 2 ^( W „) 



A 2 (x) 



X 2 



(12) 



as shown in Eq. (18) of Ref. (0) or Eq. (2.5) of Ref. 
JT3) with Aq 2 = 327r 2 e 2 iV(0)DTc, where T c is the tran- 
sition temperature of TS. In Fig. 3, the averaged value 
of A 2 , (XL = L/Io yfi( x ) ) i s plotted as a function of 
temperatures (T), where Z — 1, Rd/Rb = 1, = 1, 
Rd/Rv = 1 and Exh = 0.25A . As is shown in Fig. 3 
by curves c and d, X av > for singlet junctions. Thus 
X av is a real number and a magnetic field is screened by 
the usual Meissner effect in the DN. On the other hand, 
in triplet junctions we find A 2 „ < as shown by curves 
a and b. Therefore, the X a v becomes a purely imaginary 
number for p x + ip y - and p^-wave junctions. Formally, 
this is the consequence of the fact that the pair ampli- 
tude f n (icon , x) is purely imaginary. It is a novel feature 
of the anomalous proximity effect specific to triplet junc- 
tions that the applied magnetic field is not screened in 
the DN region and is screened only by the TS part. 
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In the lower panels of Fig. 2, the e dependence of the 
pair amplitude fx(e,L/2) is plotted for Rd/Rb = 1- For 
triplet junctions in (a) and (b), the real and imaginary 
parts of /at(£, L/2) are odd and even functions of e, re- 
spectively. On the other hand, in s-wave junctions, the 
real part of fjsr(£, L/2) is even function of e and the imag- 
inary part is odd function of e. These features are natu- 
rally understood from the fact that Real[/jv(e, L/2)] and 
Imag[/jv(£, L/2)] are respectively given by sin 9 r cosh 9i 
and cos# r sinh#i, where 6 r (9i) is the real (imaginary) 
part of 6(e,L/2). 

The unusual e dependence of /jv(e, x) in DN/TS junc- 
tions is a source of an anomalous response of DN to exter- 
nal magnetic field. When a magnetic field is applied par- 
allel to the interface, the screening current j(x) is given 
byi£ 

j(x) = ne 

N(0)DT^2 Trace[f 3 -RArOr)[f 3 , R N (x)]]A{x) 

(11) 

where A(x) is the vector potential, iV(0) is the den- 
sity of states in the normal state at the Fermi level 
in DN, uj n = 7rT(2n + 1) and T is the temperature. 



FIG. 3: Averaged value of the Meissner screening length 
X 2 av is plotted for various superconducting state with Z = 1, 
Rd/Rb = 1, Z' = 1, Rd/Ry = 1 and E Th = 0.25A , 
respectively. (&)p x + ip y -w&ve, (b)p a; -wave, (c) s-wave and 
(d)d x 2_ y 2 + id xy -wa,ve pairing. 

Let us now discuss the proximity effect in TS/DN/TS 
junctions by replacing normal electrode (see Fig. 1) by 
triplet superconductor (TS). We study the effect of an 
external phase difference tp of the p x + ip y -wave pair po- 
tentials on both sides of the junction on the LDOS in 
the DN layer. We assume that the barriers at x = and 
x = ~L are the symmetric (i.e., Rb> = Rb and Z' = Z). 
We consider the two situations, tp = and tp = ir. The 
boundary condition at x = is 

where F' 2 is given by 

4T![(1 +T?)(A' - B'fji) + 2T^B' - ^/{)] 
2 [(l+T 2 ) + 2T 1 S'] 2 -(l-T 2 )/{ 2 

(14) 
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with A' = (cos 6*0/3 _ 
cos 9 g). For ip = 0, 
and for <p = 7T, /{ : 
A o cos0/\/A2 



sin #0.9) and _B' = (sin 6*0/3 + 
we find f[ = -/ x , /£ = / 2 
= /i, /2 = -/a, with /1 = 
/ 2 



£ 2 , and ji = A sin (f>/y/ A 2 , — e 2 . As 
in the case of N/DN/TS junctions, to achieve the con- 
vergence of LDOS in numerical calculation, we replace e 
with e + 17 with small parameter 7 = 0.005Ao. 

Here, we concentrate on the LDOS at the center of 
DN, p(e, L/2), plotted in Fig. 4 for <p — and ip = tt. As 
shown in Fig. 4(a), for ip = 0, p(e,L/2) = 1 in the whole 
e range. On the other hand, for ip — tt, p(e,L/2) has a 
clear ZEP. To understand this difference, we look at the 
property of 0. By comparing the left and right bound- 
ary conditions using Eq. Q and l|14fl . we can show for 



dx 



d0_ 

dx 



Then, 



<p = that 9 L = -6» and ^ \ x =L-= ^ \x=o+- 
the resulting pair amplitudes penetrating from the left 
side superconductor and the right side interfere destruc- 
tively. Triplet pair potentials with p x + ip y -wave symme- 
try change sign under spatial inversion. The phase differ- 
ence of tt between the two pair amplitudes with MARS's 
penetrating from the two superconductors causes a de- 
structive interference in DN in the absence of an exter- 
nal phase difference. On the other hand, for ip = tt, the 
resulting 9 satisfies L = 9 and |f U=l_= -|§ |x=o+, 
and the interference becomes constructive in DN. The 
corresponding plot is also shown in the s-wave case. In 
this case, the constructive interference is enhanced for 
ip = and the resulting p(e, L/2) has a minigap. On the 
other hand, for tp — tt, due to the destructive interfer- 
ence, the resulting p(e,L/2) is unity due to the absence 
of the proximity effect at x — L/2. 



(a) p,+ip,Ap,+ipr) 



(b) s/s 




FIG. 4: Spatial dependence of p(e, L/2) is plotted for (a)p x + 
iPv/P* +iPy junctions and (b)s/s junctions with Z = Z = 1, 
Rd/Rb = Rd/Ry = 1, an d Eth = 0.25Ao, respectively for a: 
ip — and b: ip — tt. 



In this paper we have studied anomalous features of 
the proximity effect in triplet superconductor junctions 
in the framework of the quasiclassical Green's function 
formalism. We have shown that the pair amplitudes 
/iv(e) obey the unusual relation (i.e., /at(e) = —f]y(~ £ ) 
in Eq.@) in contrast to that in singlet superconductor 
junctions. This unusual energy dependence induces the 
ZEP of LDOS in DN and is also responsible for the un- 
usual behavior of a magnetic field in DN. The magnetic 
field is not screened in triplet junctions and oscillates spa- 
tially in DN. In N/DN/TS junctions, the magnitude of 
the ZEP increases with the increase of the interface resis- 
tance between N and DN which is a result of the confine- 
ment of the MARS in the DN. Therefore, to observe the 
predicted ZEP in LDOS in experiments on Sr2Ru04 and 
Ru micro inclusion systems a large resistance between 
N and DN is necessary. In TS/DN/TS junctions, the 
proximity effect is suppressed when the external phase 
difference across junctions ip = because two MARS's 
penetrating from the two superconductors cancel each 
other. On the other hand, for ip = tt, the constructive 
interference bridges the two MARS's in DN. These re- 
sults imply the phase-sensitivity of MARS. Experimen- 
tal observation of the predicted unusual proximity effect 
in triplet superconductor proximity systems would be of 
high interest. 
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